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The accumulation of defects, and in particular He bubbles, can have significant implications for the 
performance of materials exposed to the plasma in magnetic-confinement nuclear fusion reactors. Some of 
the most promising candidates for deployment into such environments are nanocrystalline materials as the 
engineering of grain boundary density offers the possibility of tailoring their radiation resistance properties. 
In order to investigate the microstructural evolution of ultrafine- and nanocrystalline-grained tungsten 
under conditions similar to those in a reactor, a transmission electron microscopy study with in situ 2 keV 
He + ion irradiation at 950 °C has been completed. A dynamic and complex evolution in the microstructure 
was observed including the formation of defect clusters, dislocations and bubbles. Nanocrystalline grains 
with dimensions less than around 60 nm demonstrated lower bubble density and greater bubble size than 
larger nanocrystalline (60-100 nm) and ultrafine (100-500 nm) grains. In grains over 100 nm, uniform 
distributions of bubbles and defects were formed. At higher fluences, large faceted bubbles were observed on 
the grain boundaries, especially on those of nanocrystalline grains, indicating the important role grain 
boundaries can play in trapping He and thus in giving rise to the enhanced radiation tolerance of 
nanocrystalline materials. 

The performance of materials in extreme environments poses important fundamental questions about the 
behaviour of condensed matter under far-from-equilibrium conditions. These conditions create challenges 
in materials design and synthesis as highlighted in a recent report 1 . Extensive research has focused on ion 
irradiation of metals 23 , semiconductors 4 and soft materials 5 to elucidate the response of these materials to extreme 
conditions. For example, the conditions found in nuclear fusion reactors have triggered recent research into the 
radiation tolerance of plasma facing components 6 , resistance to morphological changes 7 " 9 and degradation of 
their mechanical properties 10 " 12 . Athough extensive work 6,13 has focused on identifying candidate materials such 
as refractory metals (for example, tungsten and molybdenum), fundamental understanding of the atomistic 
processes which give rise to their radiation resistance has been limited by the small number of studies 14 " 16 , in 
which the dynamic response of these materials is observed directly. 

In the context of plasmas in nuclear fusion reactors, tungsten is considered one of the best options as a plasma 
facing component (PFC) material 17,18 . However, it is known that irradiation of tungsten with He (one of the main 
products of nuclear fusion reactions) can lead to microstructural changes such as bubbles 19 , pores 20 , nano- 
structures 9 and fuzz formation 21 in addition to radiation induced hardening 7 and embrittlement 13 . In the quest 
for so-called "radiation resistant" materials, efforts have included investigations into novel ultrafine- and nano- 
crystalline-grained materials 22,23 (such as tungsten in the current study) due to the conjecture that radiation 
resistance can be improved by increasing the grain-boundary area 24,25 . In the context of the current study, we 
define a radiation resistant material as one which will demonstrate reduced surface nanostructure formation 
under plasma facing conditions. Grain boundaries are known to be good He and defect sinks 26 " 28 , both of which 
can drive the aforementioned microstructural changes 29 . Research studies on nanocrystalline grained metals 
(such as copper 30 , gold 31 and nickel 30 ) and ceramics (such as Zr0 2 ) 32 have demonstrated higher radiation 
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resistance with respect to lower defect densities when compared to 
their commercially available larger-grained forms. 

Recent molecular dynamic simulations by Bai et al. 33 demon- 
strated the role of grain boundaries in annihilating interstitial defects 
and vacancies. In their work interstitials, which are highly 
mobile 15,34,35 , were shown to rapidly be absorbed by grain boundaries 
and subsequently be emitted to combine with vacancies approaching 
the grain boundary. Furthermore, enhancement of vacancy hopping 
(i.e. migration) at grain boundaries loaded with interstitials was 
observed. Grain boundaries can also trap He atoms 36 and thus 
decrease the concentration of He intragranularly. Our hypothesis 
is that an increased grain-boundary density leads to an increased 
fluence threshold (or by definition a more "radiation-tolerant" 
material) in the matrix for bubble growth, bubble coalescence and 
resulting nanostructured morphology at the surface if their forma- 
tion depends on the same as proposed in the literature 9,37 . 

To understand the effect of grain-boundary density on the radi- 
ation tolerance of tungsten and their role in trapping He, transmis- 
sion electron microscopy (TEM) investigations are invaluable 
especially when performed with in situ ion irradiation. The current 
study presents the first TEM results of in situ low-energy He + ion 
irradiation of ultrafme-grained and nanocrystalline-grained tung- 
sten. The coexistence of ultrafine (<500 nm) and nanocrystalline 
(<100 nm) grains 23 adjacent to each other enabled observation of 
the radiation response of both types of grains compared under 
exactly the same temperature and irradiation conditions. We define 
grain size by considering the distance an irradiation- induced defect 
within a grain would have to migrate in order reach its nearest grain 
boundary. For example, in a nanocrystalline grain with a nominal 
characteristic length of 100 nm an irradiation-induced defect is 
always <50 nm from a grain boundary. 

By directly observing the behaviour, size, geometry and distri- 
bution of bubbles as the He + ion fluence was increased, it was pos- 
sible to explore the role of grain boundaries in mitigating bubble 
formation in the matrix. Microstructural and morphological changes 
as observed in the current study can be correlated with enhanced 
radiation resistance of bulk ultrafine tungsten materials. This is of 
vital importance to the scientific community working on these mate- 
rials for extreme irradiating environments. This work also provides 
important observations that complement recent findings of He accu- 
mulation at semi coherent interfaces 38 . In these studies, large He 
storage concentrations are attained by the stability of He platelets 
under irradiation. Although the materials studied were different 
from our studies, the role semi coherent interfaces have on large 
He accumulation can be further elucidated by our findings. 

Results 

In order to gain a better understanding of the fundamental mechan- 
isms which govern the behaviour of defect dynamics in tungsten 
under reactor conditions, consider the 950°C and 2 keV He + ion 
irradiation conditions used in the current study. Vacancy migration 
(with an activation energy of 1.7 eV) 35 in tungsten occurs at tem- 
peratures over 800 °C whilst rapid interstitial migration (with an 
activation energy of 0.054 eV) 34 occurs at room temperature. 
Under the irradiation conditions reported here, vacancy migration, 
bubble migration and bubble coalescence can occur 29 . The particular 
novelty of this work is the investigation of grain boundary density 
engineering to control defects under extreme and far-from-equilib- 
rium conditions. The grain boundary areal density in ultrafine- and 
nanocrystalline-grained tungsten is far greater than that found in the 
microstructures of commercially available tungsten which typically 
has grain sizes that range from 1-10 um. He ions are trapped in grain 
boundaries (2D trapping) and dislocation cores (ID trapping) where 
He ions have limited diffusion compared to in the matrix (3D trap- 
ping) 27 . A large energy barrier for He diffusion back into the matrix 
makes these grain boundaries act as He sinks and effective barriers 



for intergranular diffusion of He particles. It is expected that bubbles 
will nucleate on the boundaries of grains with dimensions approach- 
ing the mean free path of migrating He and He-induced defects in 
preference to nucleating within those grains. 

Figures l(a)-(d) show in situ TEM micrographs of a sample con- 
taining both ultrafine- and nanocrystalline-grains. Grain 1 is defined 
as a nanocrystalline grain with a nominal characteristic length 
<100 nm and grains 2 and 3 as ultrafine grains with nominal char- 
acteristic lengths > 100 nm. The evolution of defects induced by low 
energy He + ion irradiation shown in Figure 1 occurred under con- 
tinuous irradiation over a fluence range of 0 to 3.2 X 10 19 ions.m" 2 . 
The TEM micrographs shown were all taken in underfocus condi- 
tions and so consequently the He bubbles appear with bright contrast 
on the inside of the bubble and thus small bubbles appear as bright 
spots. The mean number of atomic displacements generated per 
incident 2 keV He + ion was calculated using SRIM to be between 
6 and 7 events. The maximum range was calculated to have been 
about 30 nm and so within this distance there will have been suf- 
ficient numbers of vacancies for bubble formation to occur via He 
trapping by vacancies and vacancy clusters. 

Bubbles around 5 nm in size were observed at a fluence of 8 X 
10 18 ions.m" 2 . At 2.4 X 10 19 ions.m" 2 , defect clusters began to evolve 
into discernible dislocation loops in the ultrafine grains (see 
Figure 1(c) grains 2 and 3). However, in the nanocrystalline grains 
the main change was an increase in the areal density of bubbles with a 
significantly lower density of defect clusters and dislocations (see 
Figure 1(c) grain 1). Up to the end fluence of 3.2 X 10 19 ions.m" 2 , 
the mean bubble size increased in all grains (see Figure 1(d) grain 1) 
but was largest in the nanocrystalline grains. In ultrafine grains 2 and 
3, large defects clusters and dislocations were also observed (see 
Figure 1(d)) at the same fluence. 




Figure 1 | TEM micrographs of in situ 2 keV He + ion irradiation of 
tungsten at 950°C showing: (a) nanocrystalline (1) and ultrafine (2 and 3) 
grains before irradiation; (b) at a fluence of 8 X 10 18 ions.m" 2 and after 
bubble nucleation (bubbles indicated by yellow arrows); (c) after 
irradiation to a fluence of 2.4 X 10 19 ions.m" 2 showing point defect cluster 
formation (indicated by red arrows) occurred predominantly in grains 2 
and 3; and (d) after irradiation to a fluence of 3.2 X 10 19 ions.m" 2 with a 
higher areal density of point defect clusters and small dislocation loops 
evident in grains 2 and 3 whilst grain 1 demonstrates a uniform 
distribution of bubbles and a significantly lower areal density of defect 
clusters and dislocation loops, (arrows guide the eye to aid in identifying 
respective defects). 
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Figure 2 | Ultrafine and nano crystalline grains in tungsten after 
irradiation with 2 keV He + to a fluence of 3.6 X 10 19 ions.m -2 at 950°C: 
(a) overview of typical microstructure of sample; (b) higher magnification 
micrograph showing grain boundaries decorated with bubbles; 
and (c)-(d) nano crystalline grains (denoted NC) demonstrating 
significantly lower areal densities of bubbles compared to ultrafine 
grains. 

Figure 2(a) shows another selection of grains at a fluence of 3.6 X 
10 19 ions.m" 2 . The largest bubbles were observed on grain bound- 
aries (see Figure 2(b)) and in many cases appeared faceted. This can 
occur in bubbles when the size becomes sufficiently great such that it 
becomes energetically favourable to reduce the surface energy of the 
cavity by diverging from a spherical geometry 39 . The larger sizes and 
areal densities of He bubbles on grain boundaries demonstrate how 
grain size can be used as tools to engineer reduced levels of He bubble 
formation in the matrix. This will mitigate bubble coalescence in the 
matrix and/or subsequent intersection with a free surface and thus 
allow possible control over surface nanostructure formation. 
However, enhanced grain boundary grooving 40 and degradation of 
mechanical properties (for example, reduced creep resistance) 41 are 
likely to occur due to the production of large bubbles on the grain 
boundaries. This will be an important consideration in the applica- 
tion of grain size engineering to control the radiation tolerance of 
tungsten in the context of plasma facing components for magnetic 
nuclear fusion applications. 

Discussion 

To illustrate the effects of grain size on bubble formation in the 
matrix, the bubble size distributions for grains of different dimen- 
sions are shown in Figure 3. TEM micrographs showing bubble 
densities in different grain sizes are shown in Figure 4. To perform 
the bubble size distribution analysis, TEM micrographs were cap- 
tured from multiple regions and measurements made on the observ- 
able bubbles. The grain categories were determined by the grain's 
smallest dimension across its "width" and not its "length" (note: 
most grains were an elongated shape). Three grains were analysed 
for the 40-60 nm range, three grains were analyzed for the 60- 
100 nm range, and seven grains were analysed for the 100+ nm 
range. For bubble distribution analysis, the volumes of 446 bubbles 
(70 (40-60 nm range), 101 (60-100 nm range), and 275 (100 + 
range) were calculated assuming spherical geometry. These bubbles 
were randomly selected from the grains in each size category. Since 
the number of grains for each category available was low, error in 
bubble counts is given by the square root of the counts. This is a 
reasonable assumption because our counts are high and based on 
TEM observations these grains are good representations of their 
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Figure 3 | Size distributions of bubbles observed in grains of various 
dimensions at a fluence of 3.2 X 10 19 ions.m -2 . Bubbles located on grain 
boundaries were not counted, e.g. only intragranular bubbles are counted. 
The number of bubbles is counted for each grain size category and a 
standard deviation is calculated based on the number of grains investigated 
in each grain size category. The small sample size resulted in large error bars 
in some bubble size categories. 

respective grain size category. By sampling all adjacent grains in a 
region, any bias from preferential selection of grains by the observer 
were minimized. 

The bubble size distributions, which are summarized in Figure 3, 
revealed a tendency for nanocrystalline grains of less than around 
60 nm demonstrated a lower areal density compared to larger nano- 
crystalline grains (60-100 nm) and ultrafine grains (100-300 nm) 
(Figure 4). This can be understood by considering the nucleation of 
bubbles in small grains: as grain boundaries are vacancy sinks they 
reduce the probability that a migrating He atom will form a He- 
vacancy complex nucleating a bubble in the matrix; conversely, the 
probability that the He atom will instead encounter and be absorbed 
by an existing bubble will be increased. Thus smaller grains will 
exhibit fewer but larger bubbles. Figure 3 demonstrates that nano- 
crystalline grains have higher percentage of large bubbles. 
Furthermore, the grain boundaries of nanocrystalline grains with 
characteristic lengths of <60 nm were decorated with larger (with 
respect to the average bubble observed) mostly- faceted bubbles (see 
Figures 2(b)-(d) and Figure 4) than the ultrafine grains. This pre- 
dominance of larger faceted bubbles here is further evidence of 
enhanced migration (compared to commercial grain tungsten) of 
He atoms and vacancies to the grain boundaries of the nanocrystal- 
line grains. 

In addition to bubble nucleation and growth, defect cluster forma- 
tion and evolution was also observed under ion irradiation. The 
defect clusters are interstitial in nature as interstitials are in surplus 
relative to vacancies; the latter having agglomerated together with He 
atoms to form bubbles. The high binding energy of a large vacancy- 
He complex renders it stable 42 and can also lead to defect clustering 
through the process of trap mutation 43 or loop punching 8 . Trap 
mutation is defined as the process in which over-pressurized bubbles 
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Figure 4 | TEM micrographs of different grain sizes showing different bubble densities at a fluence of 3.2 X 10 19 ions.m -2 . Scale marker applies to all 
micrographs. 



emit interstitials in order to increase their volume and loop punching 
is a similar process in which an interstitial-type dislocation loop is 
generated. 

In contrast to nanocrystalline grains, some ultrafine grains showed 
a high areal density of defects (see Figure 5a). Other ultrafine grains 
featured homogenous distributions of bubbles with diffraction con- 
trast associated with them and a relatively lower areal density of 
defects (see Figure 5b). Although over-pressurized bubbles have 
never been conclusively observed under irradiation conditions, 
which induce atomic displacements 44 , we find these bubbles are con- 
sistent with loop punching and/or trap mutation mechanisms, which 
require that the bubbles be over-pressurized. 

The differences observed in the radiation response of various 
ultrafine grains are likely to be indicative of a dependence on the 
grain orientation relative to the direction of the ion beam. This can be 
due to ion beam channelling effects 45 . In a separate work, we have 
identified that grains of 001 orientation result in greater surface 
nanostructure formation than those of 110 or 111 orientation. 
Figure 5(c) shows a secondary electron image, formed using a Ga + 
ion beam, of a sample which has undergone irradiation enhanced 
recrystallization 38 ' 46 ' 47 with 30 eV He + ions at 1200°C to a fluence of 
10 26 ions.m" 2 . The faceted bubbles observed on grain boundaries 
and in grains (see Figure 2(b)) may be correlated with structures 
which form preferentially on grains of particular orientations as 
observed in experiments on bulk material irradiated at high fluence 



and temperatures (see Figure 5(c)). Figure 5(d) shows that the grain 
boundaries were bowed (ruptured upward), an observation, which is 
correlated with the high density of large faceted bubbles observed in 
the grain boundaries in this study. Further work is underway to 
investigate grain orientation effects, the correlation between the 
shape and orientation of bubbles and surface nanostructure develop- 
ment prior to fuzz formation. 

The microstructural evolution of an ultrafine grain including 
defect formation and development is shown in Figure 6 and corre- 
sponds to Grain 2 in Figure 1(a). Defect migration (see Figures 6(a)- 

(c) ), dislocation loop formation, and dislocation pinning (see Figures 

(d) and (e)) were observed. Another example of dislocation loop 
shuttling during the irradiation process is shown in Figure 7. 
Videos SI and S2 are available online in supplementary material. 
Interstitial-type defects can be trapped around He-vacancy com- 
plexes which can stabilize such high densities of defects 14 ' 48 ' 49 . The 
defects continued to evolve for about one minute after ion irradiation 
was stopped as the system returned to equilibrium (Video S3 is 
available online in supplementary material). The mobility of defects 
will have been increased by the irradiation temperature and radia- 
tion enhanced diffusion 50 effects such as the elevated point defect 
population. 

Such a high density of defects and their mobility, which were 
neglected in the nanostructure formation mechanism suggested by 
Kajita et al. 9 ' 37 may play a role in the changes that occur in tungsten in 
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Figure 5 | He + ion irradiated tungsten samples. TEM and SEM 

micrographs of ultrafine-grained tungsten taken at room temperature after 
being irradiated with 2 keV He + to a fluence of 3.2 X 10 19 ions.m -2 at 
950°C: (a) adjacent ultrafine grains with high concentrations of point 
defect clusters and dislocations; and (b) an ultrafine tungsten grain with a 
uniform areal bubble density. The diffraction contrast around the bubbles 
suggests they are combined with dislocation loops and/or are over 
pressurized as both these phenomena could induce a strain field. 
Panel (c) shows a secondary electron image formed using a Ga + ion beam 
of an ultrafine-grained sample irradiated with 30 eV He + ions at 1200°C to 
a fluence of 10 26 ions.m -2 in a complementary experiment; and (d) shows a 
cross-sectional secondary electron image of the sample in (c) showing 
nanostructured grains decorated with surface pores and a bowed grain 
boundary. 



extreme irradiating environments either directly or through their 
effects on bubble nucleation, evolution and behaviour. In a recent 
study by Ohno et al. 51 , it was demonstrated that the effect of crystal- 
lography orientation on the formation of wavy structures (ripple- 
shape structures as shown in Figure 5(c)) in tungsten exposed to 
extreme irradiating environments is related to over-pressurized bub- 
bles pushing the tungsten lattice along its slip planes. Slip can be 
enhanced by high point defect concentrations and the resulting irra- 
diation enhanced diffusion 50 which are both present under the 
experimental conditions reported in the current study. Therefore, 
the mechanisms described by Ohno et al. would be further enhanced 
by these effects. Dislocations are also trapping sites for He particles 
(ID trapping) 36 . In the current study, at the end fluence of 3.6 X 
10 19 ions.m" 2 a homogenous distribution of dislocations, on which 
bubbles were nucleated (see Figure 8), was observed in the ultrafine 
grain shown in Figure 6(f). Therefore, dislocations and dislocation 
loops can act as bubble nucleation sites and can thus affect nanos- 
tructure formation on irradiated materials. 

In summary, the application of TEM with in situ ion irradiation to 
the study of ultrafine- and nanocrystalline-grained tungsten has 
revealed several important aspects of the response of this material 
to low energy He + ion irradiation at high temperatures. The high 
radiation resistance of the nanocrystalline grains demonstrated by 
their resistance to bubble and defect formation in the matrix gives 
insight into the atomistic mechanisms, which drive these processes. 
This behaviour is correlated with the previously- reported higher 
resistance of these materials to fuzz formation which is believed to 
be initiated by bubble nucleation and coalescence. 

The larger ultrafine grains demonstrated a markedly different 
behaviour compared to nanocrystalline grains indicating the import- 
ant role of grain boundaries to the response of polycrystalline tung- 
sten to irradiation. Nanocrystalline grains demonstrated lower areal 
densities of larger bubbles compared to ultrafine grains which fea- 
tured higher areal densities of bubbles, defect clusters and disloca- 
tions. Bubbles were observed to nucleate on these dislocations 
meaning they can act as the basis of surface nanostructure formation 
thus demonstrating the importance of defect cluster and dislocation 
production on the morphology changes which can occur in He + 
irradiated tungsten. 




Figure 6 | Dynamic behavior of defects in an ultrafine grain of tungsten under continuous 2 keV He + ion irradiation at 950°C (Video SI available online 
in supplementary material): (a), (b) and (c) dislocation loop shuttling between two pinning defects (circled by solid white line) from a fluence of 
2.8 X 10 19 to 3.2 X 10 19 ions.m -2 (further example is shown in Figure 7); (d) and (e) dislocation migration and coalescence (circled by dashed white line) 
during irradiation from 3.2 X 10 19 to3.6 X 10 19 ions.m -2 ; and (f) a uniform distribution of dislocations with some decorated by bubbles at the end fluence 
of 3.6 X 10 19 ions.m -2 . 
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t = 0.000 s 
(a) 



t= 0.375 s 



3 



f = 0.500 s 
(e) 



/ = 0.625 s 
(f> 



Figure 7 | Series of consecutive video frames captured at a fluence of 2. 8 X 10 19 ions.m -2 showing a dislocation loop shuttling between two sites indicated 
by the white dashed circles in (a) and (b). The loop appears to occupy both sites in frames (d) to (f) as the shuttling frequency was greater than the video 
frame capture rate of 8 Hz. Scale marker in (a) applies to all seven panels. Video S2 is available online in supplementary material. 



The largest bubbles were observed to form on the grain bound- 
aries, rather than in the grains themselves, further highlighting the 
importance of these features as He sinks. Whilst the accumulation of 
He in grain boundaries may have significant benefits in terms of the 
reduction of surface nanostructure formation, it is important to note 
that this can also lead to degradation in the mechanical properties of 
tungsten. The engineering of the radiation resistance of tungsten via 
the control of grain size must take both these phenomena into 
account. Future work will focus on higher fluence irradiations to 
examine bubble coalescence, the fundamental process behind nanos- 
tructure formation and the effect of grain orientation. 

Methods 

Tungsten with ultrafine and nanocrystalline grains was produced using an orthogonal 
machining process which is detailed elsewhere 23 . For TEM sample preparation, small 
sections of material were mechanically polished on both sides to obtain 100 urn thick 
samples. Final thinning was performed by electropolishing with 0.1 wt-% NaOH 
solution. No variation in mass-thickness contrast, other than that expected of a 
slightly- wedged sample, was observed in the TEM across the grains suggesting that 
any preferential etching due to crystallographic orientation and/or grain size was 
negligible. This is an important consideration when comparing areal densities of 
features and phenomena, which may be affected by the proximity of free surfaces. 

Samples were He + ion irradiated in situ at 950°C within the JEOL JEM-2000FX 
TEM at the Microscope and Ion Accelerator for Materials Investigations (MIAMI) 
facility at the University of Huddersfield which is described in detail elsewhere 52 . The 
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Figure 8 | Decoration of dislocations by He bubbles: (a) dislocation 
(indicated by white circle) formed during 2 keV He + irradiation of 
ultrafine tungsten at 950°C at a fluence of 3.1 X 10 19 ions.m -2 ; and (b) 
bubbles (indicated with black arrows) observable on the dislocation at the 
slightly higher fluence of 3.4 X 10 19 ions.m -2 . 



ion beam energy was 2 keV and was operated with a flux of 3.3 X 10 16 ions.m" 2 .s _1 . 
In the MIAMI facility, the angle between the sample surface at zero tilt and the ion 
beam is 60° giving a projected range for 2 keV He + in tungsten of 10.6 nm as 
calculated by the Stopping Range of Ions in Matter (SRIM) 53 Monte Carlo computer 
code version 2013. This projected range was within the nominal characteristic length 
of the ultrafine-and nanocrystalline-grains in the tungsten samples irradiated. The 
electron energy was 200 keV and the electron beam remained incident normal to the 
sample surface during ion irradiation so that the evolution of the microstructure 
could be observed as it occurred and to allow video of the dynamic effects to be 
captured. The area monitored during the ion irradiation experiment was compared to 
adjacent regions, which were not under the electron beam (but were ion irradiated) to 
confirm there were no significant electron-beam induced effects. The sample was 
mounted and heated in a Gatan Model 652 heating holder with images and video were 
captured using a Gatan ORIUS SC200 digital camera. Imaging was performed in 
bright-field TEM conditions and images were taken with ±1.5 urn of defocus. Video 
was captured at a frame size of 480 X 480 pixels and a frame rate of 8 Hz. Additional 
post-irradiation analyses at room temperature were performed using a JEOL JEM- 
3010 TEM and an FEI-TECNAI-20 TEM with electron beam energies of 300 keV and 
200 keV respectively. 
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